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ABSTRACT 


Thermoacoustic  heat  transport  and  its  applications,  such  as  thermoacoustic 
engines  and  refrigerators,  have  been  discussed  in  a  number  of  articles  over  the 
past  several  years.  However,  lacking  from  these  articles  is  a  thorough, 
quantitative  experimental  investigation  of  the  basic  theory  underlying 
thermoacoustic  heat  transport.  Such  an  investigation  is  the  purpose  of  this  thesis. 
A  logical  starting  point  for  such  a  study  is  to  investigate  the  simplest  class  of 
thermoacoustic  engine  -  a  stack  of  short  plates  referred  to  as  a  ThermoAcoustic 
Couple,  or  TAC.  The  utility  of  this  choice  is  that  the  theory  can  be  reduced  to  its 
simplest  form  for  analysis  of  the  results.  The  results  of  measurements  of 
thermoacoustically  generated  temperature  gradients  in  TACs  subjected  to 
acoustic  standing  waves  are  reported.  The  value  of  the  temperature  gradient, 
which  results  from  an  acoustically  generated  entropy  flow  in  the  gas  in  thermal 
contact  with  the  plate,  is  a  function  of  the  acoustic  pressure  amplitude,  the  mean 
gas  pressure,  the  Prandtl  number  of  the  gas,  the  configuration  of  the  TAC,  and  its 
position  in  the  standing  wave.  Measurements  were  made  with  a  computer 
controlled  apparatus  for  drive  ratios  (the  ratio  of  the  acoustic  pressure  amplitude 
to  the  mean  pressure  of  the  gas)  from  approximately  0.1  to  2%,  in  argon  and 
helium  having  mean  pressures  from  approximately  0.1  to  0.3  MPa,  for  three 
different  TACs  as  a  function  of  their  positions  in  the  standing  wave.  The  results 
are  compared  with  predictions  based  on  a  theory  by  Wheatley  et  al  [J.  Acoust. 
Soc.  Am.  74,  153-170  (1983)].  Three  distinct  regions  of  behavior  are  apparent 
over  the  range  of  drive  ratios  investigated.  For  drive  ratios  less  than 
approximately  0.4%,  there  is  overall  good  agreement  between  theory  and 
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measurement.  For  drive  ratios  between  approximately  0.4  and  1%,  the 
agreement  diminishes  almost  linearly  with  increasing  drive  ratio.  For  drive  ratios 
greater  than  approximately  1%,  irregularities  appear  in  the  temperature 
difference  data  series  and  the  discrepancy  between  theory  and  measurement 
generally  remains  approximately  constant,  although  some  variations  occur. 
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I.  INTRODUCTION 


Thermoacoustic  heat  transport,  the  thermoacoustic  engine,  and  its 
applications  have  been  discussed  in  a  number  of  articles  over  the  past  several 
years  [Ref.  1 — 6].  Lacking  from  these  articles  is  a  thorough,  quantitative 
experimental  investigation  of  the  basic  theory  underlying  thermoacoustic  heat 
transport.  A  logical  starting  point  for  such  a  study  is  to  investigate  the  simplest 
class  of  thermoacoustic  engine  —  a  stack  of  short  plates  referred  to  as  a 
ThermoAcousticCouple,  or  TAG.  The  utility  of  this  choice  is  that  the  theory  can 
be  reduced  to  its  simplest  form  for  analysis  of  the  results.  The  purpose  of  this 
thesis  is  to  report  the  results  of  a  set  of  measurements  on  such  an  engine,  along 
with  a  quantitative  comparison  with  theory. 

The  first  measurements  of  the  thermoacoustic  effect  using  a  TAG  were 
performed  by  Wheatley,  £lal-[Ref.  1 — 3]  They  measured  the  temperature 
difference  developed  across  a  TAG  as  a  function  of  its  position  in  an  acoustic 
standing  wave.  At  relatively  low  values  of  the  drive  ratio  (the  ratio  of  the  acoustic 
pressure  amplitude  to  the  mean  pressure  of  the  gas),  the  temperature  difference 
across  the  TAG  is  a  nearly  sinusoidal  function  of  its  position  in  the  standing  wave. 
The  zeros  in  the  temperature  difference  occur  at  both  pressure  and  velocity 
nodes  (or  antinodes);  and,  the  hot  end  of  the  TAG  is  always  closer  to  the  nearest 
pressure  antinode.  Although  they  developed  a  theoretical  expression  for  the 
temperature  difference  [Ref.  6],  Wheatley,  £ial  attempted  little  quantitative 
comparison  of  their  results  with  theory.  As  they  state,  the  purpose  of  their 
measurements  was  only  to  give  them  confidence  that  they  understood  the  general 
features  of  thermoacoustic  heat  transport,  before  exploring  more  complex 
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engines.  Their  quantitative  comparisons  were  limited  to  a  single  region  in  the 
standing  wave  (near  a  pressure  node)  and  to  low  values  of  the  drive  ratio. 

All  of  Wheatley  si's  published  results  were  made  at  nearly  equal  values  of 
the  drive  ratio,  the  value  being  approximately  3  x  10*^.  However,  they  did  discuss 
the  effect  of  increasing  the  drive  ratio.  As  the  drive  ratio  "increases,  the 
[temperature  difference]  changes  from  the  nearly  perfect  sinusoid  ...  to 
something  closer  to  a  sawtooth  shape,  with  the  maxima  (minima)  shifted  to  be 
near  the  pressure  antinodes.  There  are  still  zeros  at  velocity  and  pressure 
antinodes,  but  the  slope  of  the  curve  is  very  steep  and  negative  at  the  pressure 
antinodes  [Ref.  2].  They  go  on  to  provide  a  theoretical  explanation  for  this 
behavior. 

In  light  of  this  previous  work,  we  have  undertaken  a  quantitative  investigation 
of  thermoacoustic  heat  transport.  The  results  of  the  initial  phase  of  this  work  are 
reported  here.  In  particular,  we  have  measured  the  temperature  difference 
developed  across  short  stacks  of  plates  in  acoustic  standing  waves  as  a  function 
of  the  position  of  the  plates  in  the  acoustic  field,  the  drive  ratio,  the  plate 
configuration,  the  thermal  properties  of  the  plates,  and  the  thermophysical 
properties  of  the  gas.  The  specific  aim  of  these  measurements  is  to  first  verify  the 
theory  developed  by  Wheatley  al  by  comparing  it  to  measured  results  at  low 
drive  ratios.  Secondly,  the  comparisons  will  be  extended  to  high  drive  ratios. 
These  latter  comparisons  are  especially  interesting,  because  commercially  useful 
thermoacoustic  devices  will  likely  operate  at  high  drive  ratios. 

The  thermoacoustic  effect  is  described  in  the  next  chapter,  along  with  a 
discussion  of  Wheatley's  d  al  equation  for  the  temperature  difference  developed 
across  a  TAC,  which  will  be  used  to  analyze  the  results  of  our  measurements.  A 
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description  of  the  experimental  apparatus  and  the  procedure  for  making  the 
measurements  are  explained  in  Chapter  III.  The  results  of  the  measurements  are 
analyzed  in  Chapter  IV.  The  final  chapter  contains  a  summary  of  the  thesis  and 
recommendations  for  future  work. 
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II.  THEORY 


A.  INTRODUCTION 

The  stack  of  plates  is  a  rather  simple  structure,  which  may  function  either  as 
a  prime  mover  (a  motor  producing  work),  or  as  a  heat  pump  (e.g,  a  refrigerator). 
This  chapter  is  organized  as  follows:  first,  a  description  of  the  thermoacoustic 
effect  in  a  thermally  driven  oscillator  (prime  mover)  is  given  and,  second,  the 
operation  of  a  thermoacoustic  heat  pump  is  described.  Basically,  a  large 
temperature  gradient  imposed  on  the  stack  of  plates  will  generate  an  acoustic 
standing  wave,  with  the  thermal  energy  stored  in  the  plates  being  transferred  to 
the  energy  of  the  standing  wave.  The  operation  of  the  heat  pump  is  the  functional 
reverse  of  the  prime  mover.  An  acoustic  standing  wave  imposed  on  a  stack  of 
plates  will  cause  a  temperature  gradient  to  develop,  with  the  acoustic  energy 
being  transferred  to  the  thermal  energy  of  the  plates.  After  this  discussion  an 
equation  for  the  temperature  difference  across  the  plates,  AT,  will  be  introduced 
which  was  obtained  from  Wheatley  [Ref.  2  Eq.l7]. 

B.  THERMOACOUSTIC  EFFECT 

Each  plate  is  assumed  to  have  a  sufficiently  large  temperature  gradient, VT, 
across  it  with  no  heat  conduction  in  the  longitudinal  (wave  propagation)  direction. 
Consider  a  small  parcel  of  gas  located  one  thermal  penetration  depth  from  a  plate 
and  initially  having  the  same  temperature,  T^,  as  the  portion  of  the  plate  adjacent 
to  it.  When  VT  is  large  enough,  this  gas  parcel  will  spontaneously  oscillate  at  the 
resonance  frequency  of  the  acoustic  system.  The  plates  act  as  a  source  of 


4 


energy  that  generates  the  acoustic  standing  wave.  Thus  the  plates  function  as  a 
prime  mover,  generating  acoustic  work.  The  functional  diagram  of  the  prime 
mover  is  shown  in  Fig.  1. 


output  with 
acoi'stic 
energy 


Fig.  1  Prime  mover  (r>l;  r= 


VT 


) 


The  function  of  the  prime  mover  may  be  simply  described  by  the  following 
table  : 


(1)  A  stack  of  plates  must  be  installed  in  the  resonator  and  the  material  of  those 

plates  should  be  of  low  thermal  conductivity. _ 

(2)  The  temperature  gradient  in  the  stack  should  be  sufficiently  large  and 
oriented  such  that  the  cooler  end  of  the  stack  is  adjacent  to  a  velocity 
antinode  of  the  standing  wave  and  the  hotter  end  adjacent  to  a  pressure 
antinode. 


Table  1 


Next,  the  functional  reverse  of  the  preceding  phenomena  is  considered. 

Assume  there  is  no  temperature  gradient  in  the  stack.  Acoustic  energy  is 

supplied  to  the  resonator  by  a  source  operating  at  a  frequency  corresponding  to  a 

resonant  mode.  A  simple  duct  with  a  stack  of  plates  is  illustrated  in  Fig.  2.  There 

will  be  as  many  plates  in  parallel  as  possible,  consistent  with  optimum  efficiency. 

The  plate  length  is  short  compared  with  both  the  resonant  tube  and  the  radian 

X  1 

length  of  the  acoustic  wave  X,  where  X=  ~^~ir 


Fig.  2  A  simplified  illustration  of  the  resonator  and  stack. 

The  gas  parcels  will  oscillate  with  the  driving  frequency  and  transport  heat 
from  one  end  of  the  plate  to  the  other  and  eventually  establish  a  temperature 
gradient.  A  functional  diagram  of  the  heat  pump  is  shown  in  Fig.  3. 
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Fig.  3  Heat  pump  (r<l) 


In  order  to  understand  this  effect,  the  heat  pumping  cycle  between  a  single 
plate  and  a  gas  parcel  may  be  observed  in  microscopic  view.  For  ease  of 
description,  we  will  assume  an  articulated  motion  of  the  parcel  rather  than  a 
sinusoidal  motion,  and  all  viscous  effects  are  neglected.  The  four  stages  of  the 
heat  exchanging  cycle  are  illustrated  in  Fig.  4(a)  and  4(b). 


Fig.  4(a)  Microscopic  view  of  the  heat  exchange  between  the  single  plate 

and  the  gas  particle 


V 


Fig.  4(b)  P-V  diagram  of  the  cycle  is  half  heat-exchange  &  half  without 

heat  exchange 

The  phenomena  occurring  in  a  complete  cycle  are  listed  in  Table  2. 
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Description  of  the  Phenomena 

1  -^2 

(a)  The  acoustic  energy  is  supplied  to  the  resonant  tube  and  the  gas 

parcel  is  excited  and  displaced  toward  the  pressure  antinode  by  an 
amount  Xi. 

(b)  The  compression  process  is  adiabatic. 

(c)  The  pressure  increases  from  P  to  P+AP. 

(d)  The  temperature  increases  from  Tn,  to  Tm+T* 

2-»3 

(a)  An  amount  of  heat  dQ  flows  from  the  gas  parcel  to  the  plate 

(b)  The  heat  exchange  process  is  isobaric. 

(c)  The  pressure  is  unchanged. 

(d)  The  temperature  decreases  from  T^+T'  to  Tm+T-5T 

3-^4 

(a)  The  gas  parcel  is  displaced  back  toward  the  pressure  node  by  an 

amount  Xi. 

(b)  The  expansion  process  is  adiabatic. 

(c)  The  pressure  decreases  from  P+AP  to  P. 

(d)  The  temperature  decreases  from  Tm+T-5T  to  T^-ST 

4-^1 

(a)  An  amount  of  heat  dQ  flows  from  the  plate  to  the  gas  parcel 

(b)  The  heat  exchange  process  is  isobaric. 

(c)  The  pressure  is  unchanged. 

(d)  The  temperature  increases  from  Tm-8T  to  Tm 

Table  2 


6T  is  the  temperature  change  in  the  gas  due  to  the  heat  flow  process  and  T  is 
the  temperature  change  in  gas  due  to  the  adiabatic  compression  or  expansion 
process.  The  heat  flow  dQ  is  proportional  to  the  temperature  difference  between 
the  gas  parcel  and  the  portion  of  the  plate  adjacent  to  the  parcel,  during  the  heat 
exchange  parts  of  the  cycle.  For  example,  the  temperature  of  the  plate  at  a  point 
may  be  T^+XiVT  and  the  temperature  of  the  gas  parcel  is  Tm+T.  Then  the  heat 
flux  dQ  is 
dQocAT 
and 
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AT  —  -  [Tp]ate  ■  Tg,s] 

=  -  [(T„,+XiVT)  -  (T„+T)  ] 

=  T’-XiVT 
VT 

=  T(1-— ) 

‘ ''  VT, 


cnt 


=T(i  -n 


where 


Eq.(l) 


T=-^Pi  (Ref.2Eq.4I  Eq.(2) 

and 

Tm  =  mean  absolute  temperature; 
p  =  isobaric  expansion  coefficient 
Cp  =  specific  heat  of  constant  gas  pressure 
Pm  =  density  of  gas 

Pi  =  acoustic  pressure  amplitude  (or  the  adiabatic  pressure  change) 

For  an  ideal  gas  the  product  TmP  equals  unity. 

The  value  for  VT  which  makes  AT=0  in  Eq.(l)  is  known  as  the  critical 
gradient  so  that, 

VT 

r= 

V  icnt 
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and  VTcritis  the  maximum  temperature  gradient  that  can  be  acoustically  induced 
on  to  this  plate  i.e. 

VTcrit  =  ^  Eq.(4) 

It  is  the  boundary  between  the  heat  pump  and  prime  mover. 

The  result  after  each  cycle  is  that  the  gas  parcel  has  transferred  a  small 
amount  of  heat  dQ  over  a  short  distance  X].  The  deposited  heat  is  picked  up  by 
an  adjacent  parcel  and  moved  still  further  during  the  subsequent  cycles.  The  heat 
is  shuttled  down  the  plate,  toward  the  pressure  antinode  until  it  reaches  the  end  of 
the  plate.  It  collects  there  until  it  is  conducted  back  through  the  plate  and  the 
surrounding  gas. 

When  the  cycle  operates  continuously,  it  will  heat  up  or  cool  down  the  tube  on 
either  side  of  the  stack.  That  is,  the  rigid  end  (pressure  antinode  and  velocity 
node)  will  be  heated  up  and  the  other  end  will  be  cooled  down.  The  stack  acts  as 
a  heat  engine  of  the  heat  pump  type.  The  requirements  for  the  refrigeration 
phenomena  are  listed  in  the  following  table: 

(1)  A  stack  of  plates  must  be  installed  in  the  resonator  and  the  material  of  each 

plate  should  be  of  low  thermal  conductivity. _ 

(2)  Initially,  VT=0  and  acoustic  energy  is  supplied  to  the  resonator. _ 

(3)  Because  of  the  acoustic  resonance,  the  temperature  gradient  (VT)  will 
increase,  so  that  one  end  of  the  resonator  will  become  colder  and  the  other 
end  will  become  hotter. 


Table  3 


The  entropy  flow  (AS)  occuring  in  a  complete  cycle  are  listed  as  follows: 


Process 

Description  of  the  Phenomena 

1  ^2 

(1)  Adiabatic  compression  process 

(2)  S,=S2 ;  AS=0 

2-^3 

(1)  Isobaric  heat  exchange  process 

-dO 

(2) AS-S3-S2- 

(3)  The  right  side  of  the  plate  gains  entropy  AS 

3-»4 

(1)  Adiabatic  expansion  process 

(2)  Si=  S2 ;  AS=0 

4-41 

(1)  Isobaric  heat  exchange  process 

+dQ 

(2) AS-Si-S4= 

(3)  The  left  side  of  the  plate  loses  entropy  AS. 

Table  4 


The  result  after  each  cycle  is  that  the  left  and  right  portions  of  the  plate, 
located  at  the  extreme  position  of  parcel  motion,  will  have  an  entropy  difference 
of  at  least  2AS.  During  a  single  cycle,  entropy  is  transferred  by  the  isobaric 
process  2— >3  and  4— >1 .  The  entropy  increment  per  unit  mass  from  left  to  right  is 
then  approximated  by  fRef.  2  Eq.Ma] 

=  p,(r-l)  Eq.(5) 

*m  *rn 

And  the  steady  entropy  flow  is 
S  =  pmVAS,  where  V  is  the  volume  velocity. 
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The  conflgurarion  in  Fig.  S  shows  a  solid  plate  with  length  AX,  width  y  and 

negligible  thickness.  The  length  AX  is  aligned  along  the  x-axis  and  an  ordinary 
acoustic  standing  wave  is  directed  along  X  in  the  gas  around  the  plate.  The  only 
way  that  the  heat  can  be  transported  along  X  is  by  the  hydrodynamic  transport  of 
entropy,  carried  by  the  oscillatory  velocity  Ui.  In  the  Y  direction,  the  heat-flux 
density  is  the  largest  at  about  a  distance  5k  from  the  surface  of  the  plate.  We  have 
defined  the  plate  width  as  y  because  a  cross-section  perpendicular  to  X  through 

the  plate  reveals  a  perimeter  0.  Then  the  cross-sectional  area  for  heat  flux  along 
X  is  OSk 

i.e.  (  y  *  5k)u,p  +  (  y  ♦  6k) boaom  =  OSk  [Ref.  6  p.l7  Fig.  9] 


Fig.  5  Thermal  penetration  depth  of  the  gas  around  the  plate  [Ref.  6  P.ll 

and  Fig.  9] 
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The  thermodynamically  active  volume  velocity  is  V  <«  IISkUi  Eq.(6) 
where  Ui  is  the  acoustic  particle  velocity. 

According  to  the  definition  of  the  entropy ; 


Eq.(7) 


we  then  have 


5Q  =  5ST„=p„VAST„«-p„nS«u,  4^J=^P,(r-l)T„ 

=  -n5KTn,PpiUi(r-l) 

[Ref.2Eq.l5]  Eq.(8) 


The  term 


TmP 

- —  determines  the  adiabatic  temperature  change  per  unit 

PmCp 


pressure  change  and  FTSk  is  the  volume  of  the  gas-plate  thermal  interaction.  The 
quantity  (F-l)  indicates  how  the  entropy  flow  depends  on  the  temperature 
gradient  and  that  when  VT  =  VTcrit ,  the  entropy  flow  is  zero. 

The  total  heat  flux  Q  is  proportional  to  0  but  independent  of  the  plate  length 
AX.  The  temperature  difference  AT  is  correlated  to  the  acoustic  pressure  and 
gas  parcel  velocity.  In  the  case  when  VT  «  VTeni ,  we  have 
AToepiUi 
and 

p,=po  cos(kx)  eJ"‘  [Ref.  6  Eq.  10] 

Ui=uo  sin(kx)  e'“’  Eq.(9) 
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where  k  is  the  wave  number,  x  is  the  displacement  from  the  closed  end  of  the  tube 
and  po ,  uo  are  the  magnitude  of  the  pressure  and  velocity  respectively. 

At  the  rigid  end  (i.e.  at  x=0),  pi  has  the  maximum  value  while  ui  has  the 
minimum  value.  They  are  90®  out  of  phase  with  respect  to  each  other.  Ignoring 
the  time  dependence  term  e^^in  pi  and  ui  and  multiplying  them  together  gives, 


PiUi=  pouocos(kx)  sin(kx)  =  pouo  5  sin(2kx) 

Eq.(10) 

.*.  AT  <«  poUo  sin(2kx) 

Eq.(ll) 

and  the  temperature  gradient  of  whole  section  of  the  TAC  length  (AX)  will  be: 

VT=^  Eq.(12) 

At  either  a  pressure  node  or  a  velocity  node  AT  will  be  zero.  The  relationship 
among  the  acoustic  pressure,  particle  velocity  and  the  temperature  difference  is 
graphically  shown  in  Fig.  6. 

Consider  the  fact  that  the  heat  flow  is  always  in  the  direction  towards  the 
nearest  pressure  andnode.  This  is  expected  since,  as  a  parcel  of  gas  is  moved  in 
the  direction  towards  a  pressure  antinode,  the  parcel  is  compressed  and  its 
resulting  increased  temperature  will  cause  heat  to  be  transferred  to  the  plate.  If 
the  parcel  is  moved  toward  a  pressure  node,  it  will  expand  and  its  resulting  lower 
temperature  will  cause  heat  to  be  removed  from  the  plate.  In  order  to  measure 
AT  on  the  plates,  a  thermopile  is  laminated  in  the  central  plate,  and  the  whole 
stack  is  called  a  "thermoacoustic  couple”  or  "TAC". 
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Fig.  6  The  amplitude  distribution  of  pi  ,ui  and  AT  in  the  resonant  tube 
and  their  phased  are  assigned  by  ±  sign. 

Aside  from  the  heat  flux,  the  work  flux,  i.e.  the  acoustic  power,  is  also 
important.  Ignoring  the  viscosity,  the  acoustic  work  w  is  given  by  the  area  of  the 
p-v  diagram  of  Fig.  4(b). 

w  =  PdV  Eq.  (13) 

The  total  acoustic  power  W  dissipated  by  the  plate  can  be  calculate  by 
integrating  over  all  the  gas  parcels  dxdydz  along  the  plate. 

W  =  2n6kAX-!|^P,\r-I)  [Ref.6Eq.37I  Eq.(I4) 
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So  the  acoustic  power  is  proportional  to  the  plate  length  AX  and  the  total 
energy  flow  H  is  then  the  sum  of  work  flow  and  heat  flow. 

H  =  Q  +  W  Eq.(15) 

At  the  rigid  end,  the  pressure  is  a  maximum  (called  pressure  antinode)  and  at 
the  open  end  where  there  is  a  pressure  release,  the  pressure  is  a  minimum  (called 
a  pressure  node).  A  common  term  (F-l)  exists  in  both  equations  for  Q  and  W 
given  in  Eq.(7)  and  (14). 

When  r<l  (or  when  VT<VTcrit).  then  Q  is  positive  and  the  heat  flux  is 
directed  toward  the  pressure  antinode  so  the  open  end  will  become  cooler  and  the 
closed  end  will  become  hotter.  In  this  case,  the  acoustic  power  W  is  dissipated  by 
the  gas  parcels  near  the  plates.  The  plates  function  as  a  heat  pump.  This  is 
illustrated  in  Fig.  7,  for  which 


AT=T2-Ti  =Tho,  -Teoid  Eq.(16) 

is  obviously  a  positive  value.  This  may  be  proved  by  the  measurement  of  AT.  A 
positive  half  cycle  appears  when  the  TAC  is  moving  away  from  a  pressure 
antinode,  as  shown  in  Fig.  20. 

The  prime  mover  functions  as  follows:  When  !>!  (or  when  VT>VTcrit),  then 
Q  is  negative  and  the  heat  flux  is  directed  toward  the  pressure  node.  This  tends  to 
make  the  hot  closed  end  less  hot  and  the  cold  open  end  less  cold.  In  this  case, 
acoustic  power  W  is  generated  by  the  gas  parcels  near  the  plates.  The  plates 
function  as  a  prime  mover.  This  is  illustrated  in  Fig.  8,  for  which 
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AT=  T2  -  Ti  =  Tcoid  -  Thot 


Eq.(17) 


is  obviously  a  negative  value. 


this  end 
cools  down 


Ti 


stack  of  plates 

- 


resonant 

tube 


□ 


IC 

□2 


this  end 
heats  up 


opened  end 


rigid  end 


Fig.  7  When  r<l,  the  TAC  is  a  heat  pump. 


stack  of  plates  resonant 


opened  end  rigid  end 

Fig.  8  when  r>l,  the  TAC  is  a  prime  mover 
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When  r=l  (or  when  VT=VTCTit),  no  heat  transfer  takes  place.  The  gas 
parcel  is  compressed  and  expanded  causing  its  temperature  to  increase  and 
decrease  as  it  moves.  With  this  temperature  gradient  in  the  plate,  there  is  no 
temperature  difference  between  the  parcel  and  the  portion  of  the  plate  adjacent 
to  it. 

In  our  investigation,  there  are  no  temperature  gradients  in  the  beginning. 
When  r=0,  the  heat  flow  initially  has  a  large  positive  value.  This  causes  the 
temperature  gradient  across  the  plate  to  increase  toward  the  value  of  VTcrit- 
However,  there  is  a  return  heat  flow  Qr  caused  by  this  gradient  and  the  thermal 
conductivity  of  the  plate  material,  which  is  given  by: 

Q,=kAVT  Eq.(18) 

So  the  net  heat  flow  of  the  whole  system  Qnet  will  be 

Qne,  =  Q-Qr  Eq.(19) 

=  n5KT;nPp,u,(l-n-KAVT 

At  steady  state,  Q  =  Qr ,  and  Qnci  =0,  so  VT  will  increase  only  until  a  steady 
state  is  achieved. 

The  equation  for  temperature  difference  across  the  TAC,  derived  by 
Wheatley  [Ref.  6Eq.l7] ,  is  : 
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AT  =  . 


ipjs.d+o'''^ 

(P„c(^)(l+o)) 


.  ,2o)x. 
sin( - ) 


(  i+pgs^d-g^) 

)cd  2 

(4(^)P„AX(y-1)(1-oT) 


/2fox.  V 

(l-cos(— -) ) 


where 

AT  =  temperature  difference  across  the  plate 
Po  =  peak  acoustic  pressure  of  the  driver 
6k;  =  thermal  penetration  depth 
o  =  Prandtl  number  of  the  gas 

Cp 

y  =  specific  heat  ratio  of  the  gas  (=  p  ) 


Eq.(20) 


AX  =  plate  length 

Tm  =  mean  ambient  temperature 

pm  =  mean  density  of  the  gas  [Ref.  7] 

c  =  sound  velocity  in  the  gas 
CO  =  angular  frequency 

X  =  TAC  center  position  measured  from  the  closed  end  of  the  tube 
K  =  thermal  conductivity 
d2  =  thickness  of  the  plate 


The  term  Kd2  will  be  a  function  of  the  properties  of  the  gas  and  the  plate.  For 
example,  if  a  G-10  fiber  glass  and  stainless  steel  laminated  plate  is  used,  then  the 
term  Kd2  will  be 
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Kd2  —  Kgasdgas  +  Kfgdfg+  Kg^du  +  Kgagdggg 

=  Kfgdfg+  Kgsdss  +  Kg„*(spacing) 


Eq.(21) 


Since  the  gas  path  on  each  side  of  the  plate  is  one  half  of  the  plate  separation, 
the  Kgggdgas  term  appears  twice.  This  means  the  combination  of  the  upper  and 
lower  side  gas  path  thickness  will  be  equal  to  a  plate  separation,  i.e. 


2  *  dgas  =  the  separation  distance  between  two  plates  Eq.(22) 


Measurements  were  made  on  a  TAC  that  was  mounted  on  a  probe  and 
moves  along  the  resonant  tube.  The  results  will  be  compared  to  the  theory  given 


III.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

A.  INTRODUCTION 

A  description  of  the  experimental  apparatus  is  given  in  this  chapter  along  with 
a  discussion  of  the  procedure  followed  during  data  acquisition.  The  discussion  of 
the  apparatus  will  be  divided  into  the  following  sections:  the  TACs;  the  TAC 
probe  and  associated  electrical  wiring;  the  resonator  tube  and  the  acoustic  driver 
housing;  the  TAC  positioning  system;  and  the  electronic  instrumentation.  The 
experimental  procedure  is  discussed  in  the  final  section. 

B.  EXPERIMENTAL  APPARATUS 

1.  Thermoacoustic  Couples  (TACs) 

Three  TACs  were  used  in  these  measurements.  A  description  of  the 
individual  TACs  are  given  in  Tables  5,  6,  and  7.  The  plates  of  TAC  #1  are  G- 10 
fiberglass,  while  those  of  TACs  #2  and  #3  are  a  lamination  of  302  stainless  steel 
and  G-10  fiberglass  epoxied  together.  As  shown  in  Fig.  9,  each  TAC  consisted 
of  a  stack  of  five  plates  separated  by  approximately  1.5  mm.  The  center  plate  of 
each  stack  is  a  lamination  of  two  plates.  A  thermopile,  consisting  of  several 
thermocouple  junctions  in  series,  is  located  between  these  two  plates.  The 
thermopiles  used  on  these  TACs  had  either  four  or  five  junctions.  The  purpose  of 
the  thermopile  is  to  provide  more  temperature  difference  measuring  sensitivity. 
In  order  to  determine  the  temperature  difference  developed  across  the  TAC,  the 
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thermopile  output  is  measured.  The  measured  voltage  is  divided  by  the  number  of 
thermocouples  comprising  the  thermopile  (N)  and  then  divided  by  the  sensitivity 
(in  V/°C)  of  the  particular  type  of  thermocouple  used.  Type  E,  chromel- 
constantan,  thermocouples  are  used  in  these  measurements.  Aside  from  the 
thermopile,  a  single  thermocouple  is  mounted  on  the  plate  next  to  the  center  plate 
to  allow  measurement  of  the  temperature  of  the  TAC. 


TAC#3 

Upper  layer 

•  AISI-302  stainless  steel 

•  Size;  6.85  mm  long  ♦  25.3  mm  wide  *  0.0889  mm  thick 

W 

•  Thermal  conductivity  ;  1 1.8  [Ref.  10] 

Lower  layer 

•  G-10  fiberglass 

•  Size:  6.85  mm  long  *  25.3  mm  wide  *  0.1016  mm  thick 

W 

•  Thermal  conductivity ;  0.48 

•  Thermopile  junction  pairs  number  N=5 

•  Thermocouple  wire  diameter  :  0.0254  mm  (0.001") 

•  Number  of  plates :  5 

•Thickness  of  the  TAC  :  8.15  mm 

•  Spacing  in  between  any  two  adjacent  plates :  1 .52  mm 

Table  5  The  specification  of  the  TAC#3 
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TAC#2 


Upper  layer  •  AISI-302  stainless  steel 

•  Size:  6.8  mm  long  ♦  25.3  mm  wide  *  0.0889  mm  thick 

W 

•  Thermal  conductivity :  11.8  [Ref.  10] 

Lower  layer  •  G-10  fiberglass 

•  Size:  6.8  nun  long  ♦  25.3  mm  wide  *  0.1778  mm  thick 

W 

•  Thermal  conductivity :  0.48 

•  Thermopile  junction  pairs  number  N=5 _ 

»  Thermocouple  wire  diameter  :  0.0762  nun  (0.003") _ 

«  Number  of  plates  :  5 _ 

•Thickness  of  the  TAG  :  7.7  mm _ 

•  Spacing  in  between  any  two  adjacent  plates  :  1 .52  nun _ 

Table  6  The  specification  of  the  TAC#2 

_ _ TACm _ 

Upper  layer  •  G-10  fiberglass 

•  Size:  10  mm  long  *  25.4  mm  wide  *  0.33  mm  thick 

W 

•  Thermal  conductivity  :  0.48 

Lower  layer  *  The  same  G-10  fiberglass  as  the  upper  layer _ 

•  Thermopile  junction  pairs  number  N=4 _ 

•  Number  of  plates  :  5 _ 

•Thickness  of  the  TAG  :  8.15  mm _ 

•  Spacing  in  between  any  two  adjacent  plates  :1.53  mm _ 

Table  7  The  specification  of  the  TAC#1 
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Fig.  9  The  structure  of  the  TAC  and  the  enlarged  view  of  the  central  plate 

lamination 


1  TAC  Probe 


The  TAC  must  be  moved  along  the  axis  of  a  resonant  tube  in  these 
measurements.  In  order  to  accomplish  this  positioning,  the  TAC  is  mounted  on  to 
the  end  of  a  hollow  1/8  inch  OD  stainless  steel  tube,  called  the  TAC  probe.  The 
mounting  bracket  is  illustrated  in  Fig.  10.  An  Endevco  Model  8550ml 
microphone  is  also  mounted  on  this  bracket  to  allow  measurement  of  the  acoustic 
pressure  experienced  by  the  TAC.  The  various  wires  associated  with  the  TAC 
and  the  microphone  pass  through  the  hollow  tube  and  exit  into  the  tail  section  of 
the  probe.  A  pressure  tight  feedthrough  connector  is  connected  to  the  probe  tube 
via  Swagelok  fittings.  The  feedthrough  allows  for  external  electrical  connections 
without  loss  of  pressure  in  the  resonator  tube  and  driver  housing.  The 
feedthrough  connector  is  illustrated  in  Figs.  11  and  12.  The  TAC  probe  passes 
through  a  pressure  tight  0-ring  connector  located  in  the  closed  end  of  the 
resonator  tube. 
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Fig.  10  The  structure  relationship  of  the  TAC,  microphone  and  probe 
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Fig.  11  The  anatomy  view  of  the  probe  tail  portion 
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Fig.  12  The  structure  of  the  probe  tail  portion 
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3.  Resonator  Tube  and  Acoustic  Driver  Housing 

The  purpose  of  these  measurements  is  to  investigate  the  temperature 
difference  developed  across  a  TAC  located  within  an  acoustic  standing  wave. 
The  standing  wave  is  set  up  with  a  JBL  Model — 2445J  compression  driver 
located  within  a  pressure  housing.  The  driver  housing  is  bolted  via  a  brass  flange 
to  the  end  of  a  1.217  m  long  3.8  cm  diameter  copper  tube,  called  the  resonator 
tube.  The  other  end  of  the  resonator  tube  is  also  flanged.  A  brass  plate  is  bolted 
to  this  flange  and  forms  the  closed  end  of  the  resonator.  This  brass  plate  contains 
the  connector  through  which  the  probe  tube  passes.  It  also  houses  an  Endevco 
Model  851  OB-5  high  intensity  pressure  transducer,  which  is  used  to  monitor  the 
acoustic  pressure  at  the  closed  end.  The  entire  length  of  the  resonator  is 
surrounded  by  a  7.6  cm  diameter  brass  tube.  Water  is  circulated,  with  a  Neslab 
Model  RTE-1 10  circulation  temperature  control  bath,  through  the  region  between 
the  two  tubes  in  order  to  maintain  a  uniform  temperature  along  the  resonator  tube. 
The  water  is  also  circulated  around  the  driver  housing  through  flexible  plastic 
tubing.  The  resonator  and  driver  housing  are  illustrated  in  Figs.  13  and  14.  The 
TAC  is  centered  in  the  resonator  with  a  spring-loaded  tripod  pictured  in  Fig.  15. 
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Fig.  14  The  anatomy  side  view  of  the  experiment  setup 
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4.  TAC  Positioning  System 

The  TAC  and  probe  are  positioned  within  the  resonator  with  a 
Compumotor  Model  M83-135  computer  controlled  stepper  motor  (25,000 
steps/rev.)  and  a  Compumotor  Model  2100  indexer.  The  shaft  of  the  stepper 
motor  is  connected  to  a  lead  screw,  which  in  turn  is  connected  to  the  probe 
through  a  traveling  bracket.  The  motor,  lead  screw,  and  resonator  are  attached 
to  an  aluminum  I-beam  to  provide  rigidity  and  alignment.  The  integrated 
positioning  system  is  depicted  in  Fig.  16. 
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5.  Electronic  Instrumentation 

An  instrument  controlling  and  data  acquisition  program  was  written  in 
QuickBASIC  and  was  run  by  an  IBM  AT  compatible  computer.  The  function 
generator,  the  multimeter  and  the  indexer  were  directly  connected  to  the 
computer  with  GPIB  cable  individually.  Through  the  program  command 
execution,  the  function  generator  supplies  the  acoustic  amplitude  and  frequency 
to  the  speaker;  the  indexer  controls  the  movement  of  the  stepper  motor;  and  the 
multimeter  collects  all  the  data  from  the  sensor  terminals.  An  oscilloscope  was 
used  to  monitor  the  source  and  output  of  the  speaker.  The  signal  analyzer  was 
used  to  analyze  the  spectrum  of  the  frequency  distribution,  through  which  a 
proper  resonance  frequency  can  be  picked  before  the  experiment  is  executed. 
The  temperature  reference  for  the  multimeter  was  set  up  with  an  isothermal 
block,  using  a  4-wire  method  to  measure  the  isothermal  block  temperature  by 
means  of  a  thermistor.  The  isothermal  block  is  a  common  temperature  reference 
base  for  the  tube  temperature  and  TAC.  Fig.  17  is  the  block  diagram  of  the 
electrical  data  acquisition  system.  The  temperature  reference  and  measurement 
circuits  are  shown  in  Fig.  18. 
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Fig.  17  The  electrical  data  acquisition  system 
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Fig.  18  The  utilization  of  the  TAC  and  temperature  sensors  in 

experiment  system 
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C.  EXPERIMENTAL  PROCEDURE 

Before  the  data  acquisition  is  started,  the  water  circulator  was  turned  on  and 
set  with  a  constant  temperature.  For  a  series  of  experiments  run  with  same  type 
of  gas  and  gas  pressure,  a  constant  environment  temperature  will  be  very  helpful 
for  the  analysis  of  the  results.  The  system  will  be  pumped  down  and  refilled  with 
either  Helium  or  Argon  several  times.  Approximate  resonance  frequencies  are 
determined  next  by  driving  the  JBL-2445J  speaker  with  the  noise  source  of  the 
signal  analyzer.  The  rigid  wall  Endevco  output  then  fed  back  into  the  signal 
analyzer.  The  displayed  signal  peaks  among  the  frequency  distribution  spectrum 
denote  the  resonance  frequencies  corresponding  to  the  resonant  modes  of  the 
duct.  The  JBL  speaker  input  source  is  now  changed  to  the  function  generator 
which  is  set  to  individual  resonance  frequencies  to  check  for  interference 
between  those  resonances  and  higher  harmonics  of  the  driver  signal. 

The  gain  on  the  source  voltage  amplifier  was  set  to  a  proper  range  according 
to  the  requirements  of  the  drive  ratio  (Pac/Pgas)- 

Automated  data  acquisition  now  starts.  The  program  initially  reads  voltages 
from  the  gas  pressure  transducer  and  directly  converts  them  to  pressure.  The 
reference  temperature  is  then  read  directly  from  the  multimeter  and  converted  to 
a  effective  reference  voltage  through  a  polynomial.  The  low  order  polynomial 
was  numerically  determined  from  NBS  Thermocouple  Table  over  the  limited 
temperature  range  of  interest.  Voltages  from  both  the  tube  thermistor  and  TAG 
are  read,  added  to  the  effective  reference  voltage,  and  converted  to  temperatures 
by  an  inverse  polynomial.  The  driver  voltage  amplitude  is  set  on  the  function 
generator  and  a  search  is  made  for  resonance  frequency  by  stepping  on  both 
sides  of  the  previously  determined  natural  frequency  and  reading  the  resultant 
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acoustic  pressure  amplitude.  The  frequency  producing  the  highest  acoustic 
pressure  amplitude  is  selected  and  set  in  the  function  generator. 

Once  the  preliminary  readings  are  finished  and  the  function  generator  set,  the 
data  run  commences  by  reading  the  TAG  output  voltage  and  converting  it  to  a 
temperature  as  discussed  above.  Then  the  TAG  temperature  difference  (AT) 
will  be  converted  from  TAG  output  voltage  as  following: 

=  Eq.(23) 

dT 


where 

V  =  TAG  output  voltage 

N  =  number  of  junctions  in  TAG  thermopile 

dV 

^  =  the  sensitivity  of  the  type  E  (Ghromel-Gonstantan  junction) 

thermocouple 

The  TAG  is  then  moved  to  the  next  position  with  7  mm  interval  and  a  wait 
time  of  45  seconds,  activated  and  then  readings  are  taken.  The  subsequent  data 
runs  were  made  at  constant  static  pressure  and  increasing  acoustic  pressure  by 
resetting  the  indexer  to  the  home  position  and  incrementing  the  driver  voltage 
gain  on  the  power  amplifier. 
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IV.  RESULTS  AND  DISUCSSION 


A  INTRODUCTION 

As  described  in  the  previous  chapter,  the  data  for  this  experiment  consisted  of 
measurements  of  the  temperature  difference  across  various  TACs  as  a  function 
of  their  position  in  the  standing  wave,  the  type  of  gas,  the  mean  pressure  of  the 
gas,  and  the  acoustic  pressure  amplitude.  The  results  of  these  measurements  are 
presented  and  discussed  in  this  chapter. 

B.  RESULTS  AND  DISCUSSION 

A  typical  experimental  result  for  a  low  value  of  the  drive  ratio  is  presented  in 
Fig.  19,  Two  curves  are  displayed  on  the  graph.  The  dashed  curve  represents 
the  magnitude  of  the  acoustic  pressure  (in  kPa)  as  measured  with  the  microphone 
mounted  next  to  the  TAC  as  a  function  of  the  distance  x  of  the  center  of  the  TAC 
from  the  closed  end  of  the  resonator.  The  value  of  the  pressure  magnitude  is  read 
using  the  right-hand  ordinate.  The  solid  curve  represents  the  measured 
temperature  difference  across  the  TAC  (in  °C)  also  plotted  as  a  function  of  x. 
The  value  of  the  temperature  difference  is  read  using  the  left-hand  ordinate.  In 
order  to  facilitate  comparison  with  theory,  however,  we  need  to  know  how  the 
temperature  difference  depends  on  the  position  of  the  TAC  relative  to  the 
standing  wave,  rather  than  on  the  TAC’s  position  relative  to  the  closed  end  of  the 
tube.  In  other  words,  we  need  the  temperature  difference  as  a  function  of  kx 
rather  than  x. 
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Ideally,  x  =  0  m  would  coincide  with  the  closed  end  of  the  tube  and,  therefore, 
a  pressure  antinode.  However,  if  there  is  an  error  in  measuring  the  initial  distance 
from  the  center  of  the  TAG  to  the  closed  end  of  the  tube  or  if  the  acoustic 
impedance  of  the  resonator  departs  from  that  of  the  ideal  case,  this  coincidence 
may  not  occur.  We  have  chosen  to  avoid  such  complications  by  assuming  that 
the  position  of  the  pressure  node  nearest  the  closed  end  should  correspond  to  kx  = 
2  •  The  assumed  position  of  the  pressure  node  is  then  x^ode  ~  ^  ~  ^  *  Using  the 

measured  value  of  the  frequency  and  the  thermodynamic  value  of  c,  x^o^j^  can  be 
calculated.  This  value  is  then  used  to  calibrate  the  x  measurements  by  forcing  the 
measured  position  of  the  pressure  node  to  equal  Xpo^g.  There  are  two  ways  to 
estimate  the  position  of  the  pressure  node.  The  first  is  to  measure  the  magnitude 
of  the  pressure  as  a  function  of  x,  as  displayed  in  Fig.  19,  and  then  locate  the 
position  of  the  minimum.  However,  due  to  electrical  interference  between  the 
microphone  and  the  TAG,  we  were  unable  to  measure  the  acoustic  pressure 
magnitude  and  the  temperature  difference  simultaneously.  Hence,  in  order  to 
acquire  both  quantities  two  data  runs  would  have  to  be  made  for  each  set  of 
acoustic  parameters.  In  many  circumstances,  two  data  runs  were  made  and  the 
position  of  the  pressure  node  determined  as  just  described.  However,  a  second 
method  proved  to  be  just  as  accurate  and  more  time  efficient. 

The  temperature  difference  is  proportional  to  the  product  of  acoustic  pressure 
and  particle  velocity.  Hence,  when  either  of  these  quantities  is  zero  so  is  the 
temperature  difference.  Therefore,  the  second  method  of  locating  the  position  of 
the  pressure  node  is  to  locate  the  position  of  the  zero  in  temperature  difference. 
Glose  examination  of  Fig.  19  shows  that  the  pressure  node  and  the  zero  in  the 
temperature  difference  are  not  exactly  coincident.  They  are  separated  by 
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approximately  4  mm.  Part  of  this  discrepancy  can  be  attributed  to  the  fact  that  the 

acoustic  center  of  the  microphone  may  not  be  located  at  the  center  of  the  TAG. 

In  this  case,  determining  of  the  position  of  the  pressure  node  with  the  microphone 

is  inaccurate.  After  all,  what  we  really  want  to  know  is  when  the  center  of  the 

TAG  is  at  the  pressure  node,  not  when  the  microphone  is  at  the  pressure  node. 

Also,  any  ambient  temperature  gradients  (nonthermoacoustic  in  origin)  would 

effect  the  position  of  the  zero  in  temperature  difference.  As  described  in  the 

previous  chapter,  the  resonator  was  surrounded  by  a  temperature  controlled, 

circulating  water  jacket  which  tended  to  minimize,  but  not  completely  eliminate, 

ambient  temperature  gradients.  In  any  case,  the  acoustic  wavelength  for  the 

situation  presented  in  Fig.  19  is  1.64  m,  so  even  a  4  mm  error  amounts  to  only 

0.2%  of  a  wavelength.  We  concluded  that  the  second  method  of  locating  the 

pressure  node  is  just  as  satisfactory  as  the  first.  In  the  discussions  to  follow,  we 

present  measurements  of  the  temperature  difference  as  a  function  of  kx  rather 

27Cf 

than  X.  The  product  kx  is  formed  by  multiplying  the  "calibrated"  x  by  k  =  — . 

Figure  20  illustrates  the  comparison  between  the  measured  temperature 
difference  (the  solid  curve)  and  the  theoretical  temperature  difference  (the 
dashed  curved)  for  a  drive  ratio  of  0.28%.  The  theoretical  temperature  difference 
is  computed  using  Eq.  (20).  (The  experimental  parameters  for  this  data  set,  as 
well  as  those  for  Figures  21,  22,  and  23,  are  listed  above  the  graph.)  The 
temperature  differences  are  plotted  as  functions  of  kx.  The  positions  of  the 
particle  velocity  and  pressure  antinodes  are  indicated  on  the  graph.  As  seen  in 
the  figure,  the  overall  comparison  is  good,  though  the  measured  temperature 
difference  is  approximately  20%  larger  than  the  predicted  result.  We  attribute 
part  of  the  discrepancy  to  an  uncertainty  in  the  thermal  conductivity  of  the 


42 


materials  used  to  construct  the  TAC.  Also,  as  seen  later,  measurements  taken 
with  TACs  of  essentially  the  same  design  give  slightly  different  results.  In 
general,  however,  the  measured  and  predicted  results  are  in  reasonably  good 
agreement  at  low  drive  ratios. 

Figure  21  illustrates  the  effect  of  increasing  the  drive  ratio  to  1.99%.  Four 
areas  for  comparison  are  obvious.  First,  the  agreement  is  quite  good  in  the 
vicinity  of  the  pressure  antinode.  The  slopes  of  the  measured  (solid  curve)  and 
predicted  (dashed  curve)  temperature  differences  are  approximately  equal. 
Secondly,  the  agreement  in  the  vicinity  of  the  velocity  antinode  is  poor.  The  slope 
of  the  measured  temperature  difference  is  much  less  than  that  of  the  predicted 
result.  Third,  the  data  series  of  the  measured  temperature  difference  is  quite 
irregular,  and  jagged  in  the  vicinity  of  the  maximum  temperature  difference. 
Finally,  the  measured  maximum  temperature  difference  is  significantly  less  and 
occurs  at  a  different  position  than  the  predicted  maximum.  In  short,  the 
agreement  between  theory  and  experiment  is  poor  at  high  drive  ratios  except 
near  the  pressure  antinode. 

Figure  22  illustrates  the  theoretical  temperature  difference  as  a  function  of  kx 
for  drive  ratios  from  0.17  to  1 .99%.  (The  individual  values  of  the  drive  ratio  have 
been  omitted,  in  order  to  reduce  clutter  on  the  graph.  The  data  presented  in  Figs. 
20  and  21  are  a  subset  of  these  data.)  This  graph  clearly  indicates  the  progression 
from  a  sinusoid  to  a  sawtooth  curve  and  the  displacement  of  the  maximum 
temperature  difference  toward  the  pressure  antinode  as  described  by  Wheatley 
Stal(Ref.  2].  The  actual  temperature  differences  are  shown  in  Fig.  23  for  the 
same  experimental  conditions.  The  measurements  do  confirm  qualitatively  the 
sinusoid-to-sawtooth  progression.  Notice  the  dashed  curve  demarcates  two 
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regions  of  behavior.  The  curves  betw^n  the  kx  axis  and  this  curve  are  smooth 
and  regular.  The  curves  on  the  other  side  of  the  dashed  curve  become  jagged 
and  irregular.  This  transition  occurs  in  all  of  the  data  sets  wt  have  taken. 

In  order  to  quantify  the  comparison  of  these  two  figures,  three  ratios  will  be 
examined.  The  first  ratio  is  that  of  the  experimental  slope  of  the  temperature 
difference  curve  in  the  vicinity  of  the  velocity  antinode  to  the  theoretical  slope  in 
the  vicinity  of  the  velocity  antinode.  The  second  ratio  is  the  ratio  of  the 
experimental  to  theoretical  slope  in  the  vicinity  of  the  pressure  antinode.  The  final 
ratio  is  that  of  the  maximum  experimental  temperature  difference  to  the 
maximum  theoretical  temperature  difference.  As  seen  in  Fig.  23,  the  temperature 
difference  reaches  a  maximum  on  both  sides  of  the  pressure  antinode,  so  two 
values  of  the  maximum  temperature  difference  ratio  can  be  computed.  These 
ratios  are  plotted  as  functions  of  the  drive  ratio  (in  %)  in  Fig,  24  for  the  data 
presented  in  Figs.  22  and  23.  The  slope  in  the  vicinity  of  the  velocity  antinode  is 
labeled  "SOVmax",  while  the  slope  in  the  vicinity  of  the  pressure  antinode  is 
labeled  "SOPmax".  The  maximum  temperature  difference  ratio  to  the  left  of  the 
pressure  antinode  is  labeled  "1st",  while  the  one  to  the  right  is  labeled  "2nd". 

All  ratios  have  approximately  the  same  vaiue  for  drive  ratios  less  than 
approximately  0.4%.  The  ratios  then  start  to  decrease  in  a  more-or-less  linear 
fashion  for  drive  ratios  between  approximately  0.4  to  1 .0%.  At  this  point  the 
quasilinear  decrease  stops.  It  should  be  pointed  out  that  the  data  for  the  drive 
ratio  of  1 .03%  correspond  to  the  dashed  curve  in  Fig.  23  which  demarcates  the 
regions  of  regular  and  irregular  behavior.  As  the  drive  ratio  increases  beyond 
approximately  1 .0%,  the  pressure  antinode  slope  ratio  increases  slightly  and  then 
levels  off  at  a  drive  ratio  of  approximately  1.5%.  In  the  drive  ratio  region  above 
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1.0%,  the  temperature  difference  ratio  more>or-less  levels  off,  whereas  the 
velocity  antinode  slope  ratio  tends  to  decrease,  though  at  a  slower  rate  than  in  the 
0.4  to  1 .0%  drive  ratio  region. 

Graphs  similar  to  that  in  Fig.  24  are  presented  for  several  data  sets  taken  with 
TAG  #3  in  Figs.  2S  through  27.  (The  type  of  gas  and  the  approximate  mean  gas 
pressure  is  given  above  the  graph.  Also,  since  the  two  maximum  temperature 
difference  ratios  exhibit  approximately  the  same  behavior,  the  second  one  has 
been  omitted  to  reduce  clutter.)  Each  of  these  figures  exhibit  the  same  general 
features  as  Fig,  24,  most  importantly,  the  initial  values  of  approximately  1 ,2,  the 
region  of  linear  decrease  for  drive  ratios  up  to  approximately  1 .0%,  and  the 
generally  leveling  off  behavior  beyond. 

Data  taken  with  TAG  #1  in  helium  are  presented  in  Fig.  28  and  29.  Again,  the 
same  general  behavior  is  evident.  One  noticeable  difference  is  in  the  initial  values 
of  the  slopes,  which  are  approximately  0.9  for  TAG  #1.  This  change  in  initial 
value  is  consistent  with  our  claim  that  much  of  the  discrepancy  between  theory 
and  measurement  in  Fig.  20  is  due  to  variability  between  individual  TAGs  and 
uncertainty  in  the  value  of  the  thermal  conductivity  of  the  TAG  materials. 

Data  taken  with  TAG  #2  are  shown  in  Figs.  30  and  3 1 .  Although  discussed  in 
the  previous  chapter,  TAG  #2  and  TAG  #3  are  much  the  same,  there  is  some 
variability  in  the  results  obtained  with  the  two  TAGs,  as  seen  by  comparing  Figs. 
24  and  30.  The  main  differences  are  that  the  wires  used  for  the  thermopiles  are 
different  diameters.  Although  these  data  do  not  extend  to  high  drive  ratios,  the 
same  general  behavior  is  observed  as  in  Figs.  24  through  29.  The  initial  values  of 
the  slopes  are  between  approximately  1.0  and  1.1. 
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Fig.  24  TAC#3  in  1  atm  Helium 
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Fig.  26  TAC#3  in  1  atm  Argon 
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TAC#1  In  Helium  at  1  atm 
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Fig.  28  TAC#1  in  1  atm  Helium 
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Fig.  29  TAC#1  in  3  atm  Helium 
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V.  SUMMARY,  CONCLUSIONS,  AND  FUTURE  WORK 
A.  SUMMARY 

The  goal  of  this  thesis  is  to  make  a  quantitative  investigation  of  the  basic 
theory  underlying  the  thermoacoustic  effect.  The  experimental  approach  has 
been  to  measure  the  temperature  difference  developed  across  a  thermoacoustic 
couple  (TAC)  as  a  function  of  the  position  of  the  plates  in  an  acoustic  standing 
wave,  the  drive  ratio,  the  thermal  properties  of  the  plates,  and  the  thermophysical 
properties  of  the  gas.  The  specific  aim  of  these  measurements  is  to  first  verify  the 
theory  developed  by  Wheatley  et  al  by  comparing  it  to  measured  results  at  low 
drive  ratios,  and,  secondly,  to  extend  the  comparisons  to  high  drive  ratios.  In 
order  to  make  these  measurements  a  computer  controlled  device  has  been 
designed  and  built. 

The  results  of  measurements  made  with  three  different  TACs  in  helium  and 
argon  for  mean  gas  pressures  between  approximately  100  and  300  kPa  are 
reported.  In  general,  results  can  be  summarized  as  follows.  At  low  values  of  the 
drive  ratio  (below  approximately  0.4%)  the  measurements  agree  well  with 
predictions.  The  data  series  of  the  temperature  difference  as  a  function  of 
position,  changes  from  a  sinusoid  to  a  sawtooth  as  predicted.  As  the  drive  ratio 
increases  from  0.4%  to  approximately  1  %,  the  agreement  between  measured  and 
predicted  results  decreases  almost  linearly  with  increasing  drive  ratio.  The  data 
series  continues  to  distort,  yet  remains  smooth.  At  drive  ratios  greater  than 
approximately  1%,  the  data  series  of  the  measured  temperature  difference 
becomes  irregular  and  the  linear  dependence  on  drive  ratio  stops.  At  the  highest 
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drive  ratios,  the  ratio  of  the  measured  results  to  the  jn-edicted  results  is  more  or 
less  independent  of  drive  ratio,  although  some  exceptions  were  found.  Finally, 
measured  values  of  the  temperature  difference  agree  very  well  in  the  vicinity  of 
the  pressure  antinode  regardless  of  drive  ratio.  On  the  other  hand,  they  don't 
agree  very  well  at  velocity  antinodes. 

B.  CONCLUSIONS 

Several  conclusion  may  be  drawn  from  this  work.  The  first  conclusion  is  diat 
the  thermoacoustic  effect  is  well  understood  at  low  values  of  the  drive  ratio.  In 
addition,  the  theory  developed  by  Wheatley  fiial  [Ref.  2]  agree  well  with 
measured  results.  The  second  conclusion  is  that  there  are  two  distinct  regions  of 
behavior  at  higher  drive  ratios.  One  region,  which  occurs  for  drive  ratios  less  that 
approximately  1%,  is  characterized  by  a  linear  decrease  in  the  agreement 
between  theory  and  experiment.  The  second  region,  which  occurs  for  drive 
ratios  greater  than  approximately  1%,  is  characterized  by  the  onset  of 
irregularities  in  the  temperature  difference  data  series.  Also  the  linear 
dependence  of  the  disagreement  disappears.  Further,  because  the  agreement  is 
always  good  in  the  vicinity  of  the  pressure  antinode  (velocity  node),  some  velocity 
dependent  effect  may  be  the  cause  of  the  poor  agreement  at  higher  drive  ratios. 
This  conclusion  is  also  supported  by  the  observation  of  the  sudden  onset  of 
irregularities  in  the  temperature  difference  data  series.  A  final  conclusion  is  that 
there  is  variability  between  TACs,  even  between  those  of  essentially  the  same 
design.  This  conclusion  indicates  that  the  experimental  apparatus  developed  for 
these  measurements  may  be  useful  in  searching  for  an  optimal  TAG  designs. 


56 


C  FUTURE  WORK 


One  area  for  future  work  is  to  develop  a  method  of  measuring  the  thermal 
conductivity  of  the  materials  used  to  construct  the  TACs  or  calibrating  the  TAG 
once  it  is  constructed.  A  calibration  may  be  accon^lished  by  imposing  a  known 
temperature  difference  across  the  TAG  and  measuring  the  output  voltage. 
Though  not  discussed  in  this  thesis,  we  have  made  a  fost  attempt  at  calibrating 
TAGs.  The  results  were  unsatisfactory  in  that  they  showed  wide  variability.  We 
intend  to  pursue  TAG  calibration  in  the  future.  Once  a  TAG  can  be  calibrated, 
we  can  investigate  the  causes  of  the  variability  in  agreement  between  experiment 
and  prediction  at  low  drive  ratios. 

Another  area  for  future  work  is  to  measure  the  acoustic  particle  velocity, 
most  likely  with  hot  wire  anemometry.  This  capability  will  allow  us  to  investigate 
the  conclusion  that  velocity  dependent  effects,  such  as  the  onset  of  turbulence, 
are  the  cause  of  the  disagreement  between  theory  and  experiment  at  high  drive 
ratios. 

Finally,  we  intend  to  use  the  experimental  apparatus  to  investigate  the 
behavior  of  thermoacoustic  engines  which  are  too  long  to  be  considered  TAGs, 
but  which  more  closely  resemble  the  type  of  engine  used  in  refrigerators  and  heat 
pumps.  It  should  be  possible  to  not  only  investigate  the  design  of  the  engine,  but  to 
investigate  the  optimum  placement  of  the  engine  within  the  sound  field. 
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